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Abstract 



Diffractive deeply inelastic scattering from a hadron is described in terms of 
diffractive quark and gluon distributions. If the transverse size of the hadronic 
state is sufficiently small, these distributions are calculable using perturbation 
theory. We present such a calculation and discuss the underlying dynamics. 
We comment on the relation between this dynamics and the pattern of scaling 
violation observed in the hard diffraction of large-size states at HERA. 
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The data Jl|Q] from HERA on diffractive deeply inelastic scattering can be inter- 
preted in terms of diffractive parton distribution functions, 

This function gives the joint probability per unit d^ to find a parton of type a in a hadron of 
type A and to find that the hadron is diffractively scattered. Specifically, the hadron appears 
in the final state having lost a fraction x P of its longitudinal momentum, with t being the 
invariant momentum transfer. The parton, measured at a scale fi, carries a fraction £ of 
the longitudinal momentum of the proton, or a fraction j3 — £/x F of the total longitudinal 
momentum transferred from the hadron. 

The experimental results for the diffractive deeply inelastic structure function -F^ lfT are 
related to the diffractive parton distributions by (|,|5| 

F^(x,Q^x P ,t)=j: fd ^f^^ F^/i^/^ . (2) 

q «/ Lb •-bp Lb V 

Here the hard scattering function F a is the same as in ordinary deeply inelastic scattering. 
Note that the contribution from the diffractive gluon distribution can be significant in this 
formula. Although F g is of order a s , while F q begins at order a°, one may expect that 
diffractive scattering is related to gluon exchange, so that ^> fffji- Note also that the 
diffractive parton distribution functions depend on a factorization scale // (which is typically 
set equal to Q). The \i dependence of / dlff follows from the \i dependence of F. That is, the 
diffractive parton distributions obey the ordinary perturbative DGLAP evolution equation. 
Since F is known and the evolution is known, the diffractive parton distributions at a starting 
scale fj,Q can be determined from the experimental results |2||| . 

The diffractive parton distributions at /io are not perturbatively calculable. Nevertheless, 
one would like to have some theoretically based insight into their behavior. To this end, no- 
tice that the problem lies with the large transverse size of the proton. The diffractive parton 
distributions for a physically small state would, in principle, be perturbatively calculable. 
Quarkonium, for instance, would do. Let us consider a slightly simpler model. Replace 
the proton by a special photon that couples only to a heavy quark of mass M 3> 1 GeV 
with a 7^ coupling. Then the diffractive gluon distribution can be represented in terms of 
cut Feynman diagrams such as the diagram shown in Fig. 1. The top part of this diagram 
represents the operator that defines the MS gluon distribution, as specified in more detail 
below. The lower part represents the coupling of the photon to the heavy QQ state, which 
couples again to the diffractively scattered photon in the final state. This is but one of many 
relevant diagrams. 

In this paper, we investigate this model for the diffractive gluon and quark distributions, 
taking the limit x P — > in which "pomeron exchange" dominates. Diagrams of order 
(as in Fig. 1) are the lowest order diagrams that make leading contributions, proportional 
to (l/x P ) 2 as x P — > 0. We find that we can take the limit x P — > inside the integrals that 
represent the Feynman diagrams, and express the sum of the contributing diagrams as a 
simple integral that we can evaluate by numerical integration. We thus obtain an answer 
for the diffractive gluon and quark distribution functions at a starting scale /io ~ M. 
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Evidently, the heavy quark state in the model vector meson is not the same as the light 
quark state in a proton. Nevertheless, since the model is based on nonabelian gauge theory, 
the qualitative features of the resulting diffractive parton distributions may give us hints 
about the real world. After giving the formulas that express the solution of the model, we 
explore these qualitative features. 

Consider diffractive deeply inelastic scattering of a hadron A. Let (p A , s a ) and (p A ', sa>) 
denote the momentum and spin of the incident and the diffracted hadron. Let us work in 
a frame in which the incident hadron is boosted along the positive light-cone direction and 
let p\ be its 'plus' momentum. The four-momentum transfer q = Pa~ Pa 1 has components 
qi 1 = (q + ,q~,q). (Note that in our notation q is not the virtual photon momentum, as is 
usual in the literature.) The diffracted hadron can be characterized by the fractional loss of 
longitudinal momentum x P = q + jp\ and the invariant momentum transfer t = (pa — Pa 1 ) 2 - 
For small x F , t is approximately given by t ~ — q 2 . 

Consider the definition of the diffractive parton distributions in terms of matrix elements 
of bilocal field operators 0]. This is the same definition || as for inclusive parton distribu- 
tions except that one requires that the final state include the diffractively scattered hadron. 
For gluons one has 

dx P dt 167r 2 2vr^2^y V 6 * ^ 

x(pA,s A \F a (0,y-,O) +v \p A/ ,s A r,X)(p AI ,s A r,X\F a (0)+\p A ,s A ) , (3) 

where F a (0, y~, 0) +u is the field strength operator modified by multiplication by an expo- 
nential of a line integral of the vector potential: 



F a {o,y-,o] 



Pexp \ig 1°° dx~ A+(0,x~,0) t^j 



F 6 (0,y-,Or . (4) 

ab 



The symbol V denotes path ordering of the exponential. The matrices t c in Eq. (H) are 
the generators of the adjoint representation of SU(3). The operator product in Eq. (H) has 
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ultraviolet divergences. It is understood that these are renormalized at the scale \i using the 
MS prescription. 

Similarly, for quarks of type j one has 

dfffi&x F ,t,ri ill [ dv - e -^ A v- v 

dx P dt 16 7T 2 4 7T 2 ^ J y 4^ 

x(pA,s A \%(0,y- \Q)\p A ,,s A r,X)^{p A ,,s A r,X\q j (Q)\p A ,s A ) , (5) 
where (jj(0, y~, 0) is given by 

qj(0,y-, 0) = Vexp(igJ^dx-A+{0,x-,O)t^j qj {0,y~,0), (6) 

with t c being the generators of the fundamental representation of SU(3). 

Consider now the case in which A is a photon coupling to heavy quarks, as discussed 
above. Let us consider small longitudinal momentum transfers, i P < 1, and let us examine 
leading-power contributions in x P . These contributions first appear in perturbation theory 
in graphs of the same order as the one shown in Fig. 1. We study the full set of graphs 
at this order. To perform this study, we use the lightcone gauge A~ = 0. A discussion of 
the gauge choice and the details of our analysis will be given elsewhere 0. Here we limit 
ourselves to reporting the main results and discussing their implications. To the leading 
power in x P we find that the matrix elements can be written in the following form 

(a = q,g): 



df^f A (P,x P ,q 2 ,M) 1 1 f d 2 s I f d 2 s> 



dx P dt 64 vr 2 2 V * J (2vr) 2 s 2 (q + s) 2 7 (2 vr) 2 s ,2 (q + s') 2 

xL(q,s,M,e,e')U a (x P ,[3,q,s,s')L(q,s',M,e,e') . (7) 

This result can be interpreted in terms of graphs with exchange of two gluons in a color- 
singlet state (such as the one depicted in Fig. 1). In A~ = gauge, only such graphs 
contribute leading power terms as x f — > 0. The denominators in the first line of Eq. (|7|) 
come from the gluon propagators. The factor U comes from the color-singlet projection 
of the two-gluon irreducible amplitude associated with the insertion of the operators in 
Eqs. The factors L describe the coupling of the two gluons to the incoming quark- 

antiquark system. 

More precisely, the functions U a have the form 

u a ( Xp , p, q , s, S ') = ^ p g ^2 p)x2 J 7S2 Tr k > q> s ') ^ k > q> s )) • ( 8 ) 

Here c a is the color factor and is given for quarks and gluons by 

c q = C 2 F N c , c g = C\ (iV 2 -l) . (9) 
The function u a can be written as 

u a = ^a(k, k) - ^ a (k, k + s) + 4> a (K k - q) - ^ a (k, k - q - s) . (10) 
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For a = q, the function if) has the following expression in terms of the 2x2 Pauli a matrices: 



U> k, p) = -o ; ; -P'ff . 11) 

For a = g, the function if) is expressed in terms of two transverse vector indices i,j = 1, 2 as 
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The functions L in Eq. (|7|) are given by 

L ^ M ^= e WlwkJ° dz (13) 

xTrj -&(z,r + s + zq, M,e') -&(z,r - s- (1 -z)q,M,e') 

+ $^2, r + z q, M, e') + $t(z, r - (1 - z) q, M, e')] $(*,r,Af,e)} , 

where eg is the quark electric charge in units of e = y/Ana and 

$(z,k,M,e) = , - 1 r- f(l - z) e ■ ak ■ a - zk ■ oe ■ a + i M e ■ a] , (14) 

(k 2 + M 2 ) 

with e and e' being the initial and final photon transverse polarizations. 

In Eq. (0), the Green functions U a are universal: as long as the diffracted system is small 
enough that lowest order perturbation theory applies, these functions control the scattering. 
The dependence on the specific diffracted system is contained in the functions L. In Ref. [0, 
we will utilize a selection of hadronic functions L. 

The explicit form of the Green function U a is given in Eqs. (^)-(jl^). In the framework 
of x~-ordered perturbation theory, the function if) in Eq. (|10D can be interpreted J?]-^| as 
the wave function for an effective quark or gluon state created by the operator in Eq. (||) or 
Eq. (0). (Here the line integral of A + in Eqs. (f§),(|6]) represents the interaction of an effective 
infinite-momentum parton with the gluon field.) The wave function description is natural in 
approaches that look at diffractive scattering in a frame in which this quark or gluon system 
has large 'minus' momentum PJTOfl. Note that each of the terms in Eq. ([lOj) would give rise 



to an ultraviolet-divergent integration over k in Eq. (|S|), but that the bad behavior cancels 
among the terms. This is because k 2 — ► oo corresponds to the partons being at the same 
transverse position. But since the net color of the state is zero, the coupling of the state to 
gluons vanishes in this limit. 

At the leading power level, the x F dependence of the diffractive parton distributions is 
given simply by the overall factor that we have factored from the integral in Eq. (^j). 
The (3 and t (or q 2 ) dependences, on the other hand, are non-trivial and come from the 
factors in the integrand in the right-hand side of Eq. (|7]). The only explicit (3 dependence 
in this integrand, in particular, is contained in the functions U a . Inspection of Eqs. (§)- 
(O) shows that, for (3 — > 0, the quark and gluon Green functions behave respectively as 



U q ~ const./? and U g ~ const./? 1 . For ft — > 1, the functions U a have a constant behavior 
(1 - ft) at finite q 2 , s 2 and s' 2 . 
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We now evaluate the diffractive distributions by performing numerically the integrations 
(0) ; © ) (HD • 111 Fig. 2 we report the results by plotting the /3-dependence of the quark and 
gluon diffractive distributions for different values of q. To emphasize the regions of small 
(3 and large (3 we make a logarithmic plot in the variable (3 /{I — (3). The curves in Fig. 2 
reflect the behavior of the functions U a discussed above. In particular, as (3 —>■ 1 the 
distributions have a constant behavior. The asymptotic constants, on the other hand, are 
small compared to the values of the distributions at intermediate (3. Correspondingly, the 
diffractive distributions fall off as one approaches the small (1 — f3) region. Note that the 
gluon distribution is much larger than the quark distribution. Roughly, the different order 
of magnitude is accounted for by the color factors in Eq. (|9|), c g /c q = 27/2. 
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FIG. 2. 
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The (3 dependence of the gluon (above) and quark (below) diffractive dis- 



tributions for different values of q 2 ~ \t\. The rescaled distributions h a are defined as 



x 2 P M 2 



df*f A /{dx P dt)\ /{a 2 e%a 



The calculation described so far does not have scaling violation. One may interpret the 
results above as a model for the diffractive distribution functions at a scale of order M 2 
and study the scale dependence due to renormalization group evolution. In Fig. 3 we report 
the result of using ordinary evolution equations for the evolution of the initial distributions 
up to different values of Q 2 . In this figure we plot the flavor-singlet quark distribution. In 
leading order this is proportional to the structure function F$ ■ Here we have assumed the 
initial scale and the mass to be Ql = M 2 = 2 GeV 2 (about the charm quark mass squared) 
and we have integrated the distributions df^/[dx F dt] over t ~ — q 2 from to M 2 . 

Is there any relation between these calculations and the data from diffractive deeply 
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FIG. 3. Q 2 -evolution of the diffractive singlet quark distribution 

£(/?, Q 2 ) = M J2j So dq 2 dfjf^/ (dx P dt) , where Af = x F /{a 2 eQO^) and the sum runs over 
j = {u,u,d,d,...}. 

inelastic scattering at HERA? Obviously, the protons probed in experiments at HERA have a 
large transverse size, in contrast to the small-size hadronic state considered in our calculation. 
Nevertheless, we find, by studying different hadronic wave functions , that the /5-behavior 
as well as the relative size of the singlet quark and gluon distributions depends only weakly 
on the specific nature of the incoming hadronic state. This suggests that qualitative features 
of the results for small-size hadrons may apply with more generality. 

In particular, one of the most peculiar features of the HERA data is the striking 
difference in the Q 2 evolution between the diffractive and the inclusive structure functions. 
In our calculation, the diffractive quark distribution grows with Q 2 at low /3 and decreases 
at high f3 (Fig. 3). The stability point at which the behavior changes is /3 ~ 0.5. This is 
in striking contrast with the case of the inclusive quark distribution in a proton, for which 
the stability point is at x ~ 0.08, but is in qualitative agreement with the behavior found 
in the HERA experiments. The explanation for growth in the quark distribution up to such 
large values of f3 in the diffractive case is that the initial gluon distribution (Fig. 2) is large 
even at large values of (3. As Q 2 increases, the gluons feed the quark distribution through 
the splitting g — > qq. 

This research is supported in part by the US Department of Energy grant DE-FG03- 
96ER40969. Z.K. acknowledges the hospitality of the Institute of Theoretical Science at the 
University of Oregon and the Theory Division at CERN while part of this work was being 
done. 
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